In [2]:
using Plots,ApproxFun,Interact,Reactive; gadfly();
Evaluate the first cell to start the plot, then the second cell to evolve the plot. Evaluating the second cell will continue evolving the solution. It may take two times before it works.
In [1]:
u0=Fun(x->2exp(-x^2),[-10.,10.])
d=domain(u0);D=Derivative(d)
g(y)=-6.y*y'
x=Signal(u0);map(y->plot(y;ylims=(-1.,3.)).o,x)
In [3]:
u0=BDF2([neumann(d);rdirichlet(d)],-D^3,g,zeros(3),u0,0.004,300,x,10E-7);